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摘 要 : 目前 电动 汽车 常 以 直流 无 刷 电 机 (BLDCMD 作 为 驱动 器 , 但 BLDCM 调 速 控制 系统 中 模糊 控制 器 的 量化 因子 和 
比例 因子 采用 传统 方法 自 调节 能 力 弱 ， 针 对 该 问题 提出 一 种 改进 QPSO 算法 (AMF-QPSO) 实 现 对 量化 因子 和 比例 因子 
的 自 适 应 调节 。AMF-QPSO 算法 以 收缩 一 扩张 系数 (contraction expansion，CE) 控 制 方式 为 研究 重点 ， 提 出 粒子 活性 概 
念 ， 并 以 其 作为 反馈 量 ， 实 现 动 态 自 适 应 调节 CE 系数 。 同 时 ， 为 防止 种 群 高 度 聚 集 ， 采 用 精英 群体 随机 交叉 学 习 机 
制 , 对 部 分 活性 低 的 精英 粒子 进行 扰动 , 增强 种 群 后 期 多 样 性 。 最 后 , 通过 LabVIEW 实验 平台 , 以 具体 案例 验证 AMF- 
QPSO 算法 性 能 。 实验 结 果 表明 ,AMF-QPSO 优化 的 模糊 PID 控制 器 具有 比 标准 模糊 PID 控制 器 和 QPSO 优化 的 模糊 
PID 控制 器 更 好 的 控制 性 和 自 适应 性 。 
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Abstract: At present, electric vehicles often use DC brushless motor (BLDCM) as the driver, but the traditional method is 
difficult to design quantization factor and proportion factor of fuzzy controller in the BLDCM speed control system. An 
improved QPSO algorithm (AMF-QPSO) is proposed for the on-line tuning of the quantization factor and proportion factor. The 
AMF-QPSO algorithm takes the contraction-expansion coefficient (CE) control method as the research focus, proposes the 
particle activity concept, and uses it as the feedback quantity to realize the dynamic adaptive adjustment CE coefficient. At the 
same time, in order to prevent the high aggregation of the population, the elite group random cross learning mechanism is used 
to disturb Some of the elite particles with low activity and enhance the diversity of the population in the later period. Finally, 
through the LabVIEW experimental platform, the performance of the AMF-QPSO algorithm is verified by a specific case. The 
experimental results show that the fuzzy PID controller optimized by AMF-QPSO has better control and adaptability than the 
standard fuzzy PID controller and the fuzzy PID controller optimized by QPSO. 
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expansion 
引言 势 ， 受 到 了 伺服 控制 、 交 通 运输 、 物 料 输送 、 工 业 自 动 化 等 领 
加 或 人 员 的 青睐 上。 在 交通 运输 领域 ， 新 能 源 技术 的 推广 以 及 绿 


近 些 年 , 随 着 新 材料 、 先 进 制造 、 电力 电子 等 技术 的 发 展 ， 色 崛 起 战略 的 实施 ,以 BLDCM 为 驱动 器 的 电动 汽车 日 益 普 及 。 


前 现 出 了 一 批 怕 


ys 


~ 


PT 


能 优越 的 新 型 电机 。 无 刷 直流 电机 (brushless 然而 电动 汽车 调 速 控 制 系统 具有 非 线性 强 、 变 量 多 、 耦 合 程度 


direct current motor, BLDCM) 作 为 新 型 电机 典型 代表 ， 由 于 损耗 高 、 控 制 器 参数 难 整定 等 特性 ， 提 高 系统 可 控 性 和 自 适 应 性 ， 


小 、 功 率 密 度 高 、 过 载 能 力 强 、 噪 声 小 、 可 控 性 强 等 一 系列 优 多 化 参数 整定 方式 已 成 为 进一步 推广 电动 汽车 的 重要 和 手段。 
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录用 稿 


模糊 PID 控制 器 作为 先进 控制 器 典型 代表 , 运用 模糊 数学 理论 ， 
将 高 等 智能 生物 模糊 语言 转换 为 机 器 可 认 知 的 程序 语言 ， 算 法 
见 机 制 类 似 于 人 类 思维 模式 ， 可 在 被 控 对 象 非 完全 可 知 条 件 
空 制 ， 适 用 于 电动 汽车 调 速 控制 系统 设计 与 参数 优化 问 


但 在 实际 工业 应 用 中 ,模糊 PID 控制 器 参数 缺乏 在 线 调整 
能 力 ， 造 成 系统 控制 精度 不 高 、 抗 干扰 不 强 。 随 着 智能 算法 的 
发 展 ， 将 智能 算法 应 用 于 模糊 控制 器 参数 整定 已 成 为 研究 热点 
5。2004 年 , Sun 等 人 四 受 量子 空间 中 粒子 启发 , 提出 了 QPSO 
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1 ”量子 行为 粒子 群 算法 


QPSO 算法 是 在 PSO 算法 基础 上 提出 的 一 种 群 智能 算法 ， 
该 算法 摆脱 传统 速度 -位 移 轨 道 模 型 限制 9, 采用 蒙特 卡 罗 算 法 


确定 粒子 更 新 后 位 置 ， 算 法 智能 化 程度 更 高 。 其 进化 方程 如 式 
(1) 所 示 。 
X(t+D)=p, zc u(t) e U(0,1) (1) 
WW u(t) ' 


其 中 : Pi 为 局 部 吸引 子 , 决定 粒子 寻 优 方向 ， 5 为 势 阱 特 


算法 。 量 子 空间 中 每 个 粒子 位 置 相当 于 一 个 可 行 解 ， 其 位 置 坐 
标 由 波 函 数 决 定 ， 且 状态 遵循 琶 加 原理 ， 具 有 很 强 随 机 性 ， 智 
能 化 程度 高 。 
标准 QPSO 算法 尽管 有 许多 优点 ,但 依然 存在 不 足 ,其 中 ， 
收缩 一 扩张 系数 (contraction-expansion，CE) & 是 除 种 群 规模 、 
粒子 维度 以 及 最 大 迭代 次 数 外 ， 唯 一 需 人 工 设 定 的 控制 参数 ， 
其 控制 策略 对 算法 性 能 影响 大 四。 文献 [8] 证 明 CE 系数 必须 满 
足 a<1.781 才 能 使 粒子 可 靠 收 敛 ， 且 比较 了 固定 系数 策略 和 非 
线性 递减 策略 优 缺 点 。 文 献 [9] 提 出 了 动态 非 线性 递减 CE 系数 
空 制 策略 ， 根 据 不 同 优化 问题 特性 采用 凹凸 性 不 同 控制 曲线 。 
但 无 论 是 固定 值 、 线 性 递减 ， 还 是 非 线 性 递减 CE 系数 控制 策 
略 ， 均 属 无 反馈 控制 策略 ， 系 数值 及 控制 曲线 由 仿真 实验 或 经 
验 公式 求 得 ,这 类 系数 控制 策略 最 大 不 足 表 现在 自 调节 能 力 弱 ， 
算法 只 能 按照 设 定好 的 调节 方式 调节 CE 系数 值 ， 未 能 根据 粒 
子 搜索 阶段 的 转变 而 调整 系数 控制 方式 。 鉴 于 此 ,研究 人 员 提 
出 了 具有 不 同 反馈 机 制 自 适应 系数 控制 策略 。 文 献 [10] 利 用 粒 


子 距 全 局 最 优点 的 距离 作为 CE 系数 调控 依据 ,文献 [11] 提 出 了 
进化 速度 因子 和 聚集 度 因子 的 概念 ， 用 于 在 线 修正 粒子 CE 系 
数值 。 


受 上 述 文献 启发 ， 本 文 在 探究 CE 系数 与 势 阱 长 度 关系 基 
础 上 ， 提 出 粒子 活性 概念 ， 以 势 阱 长 度 相对 变化 率 作 为 粒子 活 
性 度量 标准 , 并 以 粒子 活性 作为 反馈 量 , 实现 自 适 应 调节 CE 系 
数 ， 与 此 同时 ， 为 提高 收敛 精度 ， 防 止 粒子 陷入 局 部 最 优 ， 提 


征 长 度 ， 其 大 小 直接 制约 粒子 搜索 范围 ，u() 为 在 (0,1) 上 的 随 
机 分 布 。Clerc 等 人 0 对 PSO 算法 中 粒子 收敛 行为 研究 表明 ， 


要 保证 算法 收敛 ， 必 须 使 每 个 粒子 i 均 收 敛 到 各 自 局 部 吸引 子 
Pij 点 ， 其 中 吸引 子 坐 标 如 式 (2) 所 示 。 
0 oy 和 


Cj(D + en CD 


pij(D) =@, DP,(D)+ 1 一 多， GO]c (D,1<7<N 
cnij(D (3) 


VD = An D+ oon) D+ on 


其 中 : 和 < 分 别 为 粒子 局 部 学 习 因 子 和 全 局 学 习 因 子 ， 00 
入/ 为 (0 上 的 随机 分 布 ， ?0 为 粒子 ;第 了 了 维 的 局 部 最 优 
位 置 ，6,0 为 粒子 i 第 / 维 的 全 局 最 优 位 置 。 
要 保证 粒子 收敛 至 各 自 吸引 子 ， 则 必须 满足 

lmij0=0 (4) 
因此 , 如 何 控制 50 成 为 确保 算法 收敛 以 及 影响 算法 性 能 
的 关键 典型 50 控制 方式 为 利用 粒子 距 群 体 平均 最 优 位 置 的 
距离 进行 调控 ， 记 种 群 平均 最 优 位 置 坐标 为 co ， 其 坐标 可 用 
式 (5) 表 示 。 


1 之 
CO=(GO,CO…CoO)= 方 之 BO) 
Fl 
1 1 之 1 妃 
RO SRO A) 


其 中 : D 为 粒子 维度 ; WN 为 种 群 规 模 。 此 时 ， 粒 子 势 阱 长 度 可 


(5) 


出 精英 群体 随机 交叉 学 习 机 制 ， 使 得 活性 低 的 精英 粒子 有 机 会 
与 种 群 中 其 他 粒子 交换 维度 信息 ， 跳 出 局 部 最 优 ， 搜 寻 到 全 局 
范围 更 优 解 。 精 英 群 体 随机 交叉 学 习 机 制 既 保 证 了 种 群 收敛 速 
度 ， 又 在 一 定 程度 上 提高 了 种 群 后 期 多 样 性 。 最 后 ， 针 对 电动 
汽车 调 速 控制 系统 中 模糊 PID 控制 器 的 量化 因子 和 比例 因子 在 
线 调节 能 力 弱 的 问题 ， 应 用 改进 QPSO 算法 (adaptive mutation 
QPSO algorithm based on feedback mechanism, AMF-QPSO) 进 行 
优化 设计 。 为 验证 算法 性 能 ， 在 分 析 实 际 应 用 场景 基础 上 ， 通 
过 LabVIEW 实验 平台 ,以 具体 控制 案例 验证 AMF-QPSO 算法 
性 能 , 与 基本 模糊 PID 及 标准 QPSO 算法 优化 后 的 模糊 PID 控 
制 器 进行 性 能 对 比 ， 实 验 结果 表明 改进 后 的 控制 器 具有 更 好 的 
控制 性 和 自 适 应 性 。 


(6) 式 表示 。 


ne 2:.0:|C)— xX, (6) 
其 中 : & 即 为 CE 系数 。 将 式 (6) 代 入 式 (1)， 得 最 终 进 化 方程 如 
(7) 式 所 示 。 


1 
X(t+D= pM)+alc,D) -XD In 司 (7) 


2 ”基于 反馈 机 制 的 自 适 应 变异 QPSO 算法 


2.1 CE 系数 与 势 阱 长 度 之 间 / 的 关系 研究 

式 (9 可知 ， 势 阱 长 度 调节 过 程 中 CE 系数 & 起 着 决定 
性 作用 ，< 取 值 不 同 ， 最 直接 的 效果 是 增 大 或 者 削减 粒子 搜索 
空间 。 当 x 取 较 大 值 时 ， 势 阱 长 度 L 相对 较 大 ， 粒 子 具 有 较 强 
全 局 搜索 性 能 ， 但 当 x 过 大 时 ， 会 造成 种 群发 散 ， 当 “ 取 较 小 
值 ， 势 阱 长 度 L 相对 较 小 ， 粒 子 具有 较 强 局 部 搜索 性 能 ， 但 当 


过 小 时 ， 粒 子 将 失去 搜索 性 能 ， 种 群 进化 处 于 停 济 


想 的 CE 调控 方式 为 前 


带 状 态 。 理 


期 大 ， 使 粒子 具有 较 强 全 局 搜索 能 


后 期 小 ， 使 粒子 具有 较 强 


ta 


虽 局 部 搜索 能 力 。 
因此 ， 在 保证 算法 收敛 前 提 下 ， 确 定 势 B 
[2 


长 度 工 控制 策略 


关键 在 于 选择 CE 系数 


的 控制 方式 ， 


故 研究 & 对 于 势 阱 长 度 


Z 的 影响 ， 探 究 两 者 间 的 关系 对 于 算法 性 能 改进 研究 具有 很 好 


指导 作用 。 为 便于 分 析 , 对 一 维 情况 下 的 进化 公 


式 (7) 进 行 化 简 ， 


将 ? 点 和 C 均 固 定 在 原点 , 令 粒 子 初始 位 置 X(0)=1, 则 简化 后 


的 进化 公式 如 (8) 所 示 。 


XUC+D=ewlXONnCL) 
ui(D) 


由 (8) 式 可 知 , 要 保证 
则 此 时 工 需 满足 


(8) 


算法 收敛 , 即 粒 子 收敛 至 > 点 (原点 )， 


limL=limalX,()|=0 (9) 

式 (9) 可 知 , 算法 是 否 能 可 靠 收 敛 可 转换 为 观察 L 是 否 收 

敛 。 因 此 ， 根 据 式 (8) 进 行 随机 模拟 实验 。 图 1 给 出 了 x 取 不 同 
值 时 ， 势 阱 长 度 Z 的 变化 情况 。 其 中 横 坐 标 为 迭代 代数 上， 纵 坐 
标 为 势 阱 长 度 L 。 为 更 加 清晰 观察 两 者 间 关 系 , 取 以 10 为 底 的 


对 数值 表示 ， 即 log(Z) 。 


1 x10’ 


0 2.4 
1 x 10 
100 
-80 w=15 
全 -260 
名 -440 
-620 
3000 09 18 27 36 4.5 
1 x10’ 
20 
-20 
-60 
古 
2 -100 
-140 
-180 
0 0.4 08 12 1.6 2.0 


1 x 10 


1.8 2.7 3.6 4.5 


48, 2 6 12 


图 1 CE 系数 ca 与 势 | 


0 14 28 42 5.6 7.0 
X 10? 


[em] 
So 
2 
So 
00 
= 
MD 
一 
CN 
MD 
SS 


0 0.2 0.4 0.6 0.8 1.0 
1 x10! 


长 度 工 的 关系 


表 1 记录 了 cx 取 不 同 值 时 ，L 
者 工 不 收敛 时 的 最 大 办 代 次 数 。 


改 敛 时 的 最 大 迭代 次 数 了 或 


表 1 CE 系数 & 及 相应 迭代 代数 
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[24 0.3 0.6 0.9 1.2 13 | 1.79 2.0 
T 500 800 1200 2000 6000 10000 20000 20000 
由 图 1 和 表 1 可 知 ， 当 w<177 时 ， 随 种 群 更 新 ， 势 阱 长 度 


逐渐 减 小 ， 
趋势 ， 算 法 
出 现 波动 ， 
时 
势 


算法 不 收 伍 。 上 述 实验 结果 论证 


且 需 要 长 时 间 的 种 群 
, 势 阱 长 度 收敛 过 程 波动 更 为 


直至 为 零 ， 表 明 这 一 参数 范围 内 势 
色 敛 ; 当 < 小 于 但 接近 


1.78 时 ， 势 阱 


直 呈 现 收缩 
长 度 收敛 过 程 


的 高 度 关 联 性 ， 且 当 < 取 不 同 值 时 ， 势 B 


更 新 才能 收敛 至 零 ， 当 w>1.79 
剧烈 ,而且 经 20 000 次 迭代 后 ， 
阱 长 度 仍然 很 难 收敛 至 零 ， 说 明 粒 子 难以 收敛 至 P 点 ， 此 时 
了 CE 系数 与 势 阱 长 度 工 之 间 
长 度 L 的 收敛 速度 与 


收敛 精度 均 存 在 区 别 ， 间 接 造 成 算法 收敛 速度 和 收敛 精度 的 不 


同 。 


2.2 具有 反馈 机 制 的 CE 系数 自 适应 控制 策略 


l 


才能 使 算法 收敛 ， 且 对 于 
度 工 的 控制 。 


QPSO 算法 中 ， 工 决定 粒子 


由 2.1 节 仿 真 结果 可 知 ，CE 系数 只 有 在 合适 范围 


内 取 值 ， 


CE 系数 的 调控 实质 上 是 对 于 势 B 


长 


变化 率 则 反映 J 粒子 搜索 范围 


群 察 集 度 高 ， 粒 子 较 为 集中 ， 此 时 势 阱 机 


变化 


I 


搜索 范围 
情况 。 根 据 式 (1) 可 得 ， 若 种 
对 变化 率 将 很 小 或 者 


， 工 的 相对 


不 再 变化 , 说 明 粒 子 已 收敛 至 吸引 子 , 更 新 基本 处 于 停滞 状态 ; 


若 种 群 聚集 程度 低 ， 


进行 寻 优 。 因 此 ， 可 用 势 
粒子 活性 因子 ww 定义 如 


Si Os 


粒子 较为 分 散 ， 此 
化 率 较 大 ， 说 明 粒 子 位 置 更 新 处 于 
长 度 术 
式 (10) 所 示 。 


OO 一 ACE- 引 


活跃 状态 ， 


时 粒子 势 


粒子 还 在 大 范围 


长 度 相对 变 


EE) 


其 中 :5,0 为 粒子 i 第 j 维 更 


新 至 


为 粒子 :第 7 维 


更 新 至 +t-3 代 时 的 


维 更 新 至 :时 的 活性 因子 。 


势 


于 算法 收敛 时 , 势 P 


: 代 时 的 势 了 
长 度 ， 


长 度 ; 
5 为 粒子 i 第 j 


日 对 变化 率 来 表示 粒子 活性 。 


(10) 


L(t—3) 


呈现 收缩 趋势 ， 


势 阱 长 度 逐 代 递 减 ， 


L,(D) <L,(t—3), 


经 理论 分 析 可 知 ,，s 
说 明 粒 子 活性 越 强 ， 


粒子 
所 以 sj(D e[0， 
j 越 大 , 说 明 势 阱 长 度 3 代 内 收缩 剧烈 ， 


粒子 向 全 局 最 优 位 置 移动 速度 越 快 ; 


也 。 


向 全 局 最 优 位 置 移动 ， 即 


小 ， 说 明 势 
向 全 局 最 优 位 置 移 动 速度 慢 。 


2.1 节 探讨 了 CE 系数 与 势 B 
影响 的 取 值 ， 


表明 CE 系数 直 


分 析 可 得 ， 固 定 系数 、 线 性 或 
策略 都 是 一 种 无 反馈 的 系数 控 
控 量 工 之 间 的 关系 。 基 于 这 


制 的 CE 系数 控制 策 
5 作为 反馈 量 ， 构 建 
定义 如 (11) 式 所 示 。 


闭环 


长 度 工 之 


间 的 关系 ， 


5 越 


长 度 3 代 内 收缩 缓慢 ， 说 明 粒 子 活性 越 弱 ， 粒 子 


仿真 结 


制 策 


发 现 ， 
各 ， 以 工 的 相对 变化 率 ， 即 粒子 活 
空 制 模型 ， 实 现 %j 的 动态 调整 。 其 


制约 着 粒子 搜索 范围 。 而 且 


FE 线性 递减 以 及 部 分 
作 ， 并 未 考虑 CE 


适应 控制 
系数 与 被 


本 文 提 出 


种 具有 反馈 机 


Cij(D = Co +(1—s;,(D) 4 


其 中 :，w 为 CE 系数 初 值 ; 
活性 因子 


4 为 (0.0 中 随机 数 ; 
0 与 $j 成 负 相 关 。 当 sj 越 小 ， 


生 因 子 


(GD 
sj(D 为 粒子 i 


代表 粒子 搜索 能 力 


201810.00039v1 
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录用 稿 


越 弱 ， 此 时 应 适当 增 大 “,， ， 维 持 种 群 多 样 性 ; 
表 粒 子 搜索 能 力 越 强 , 此 时 ww, 适 
步 长 过 大 ， 越 过 最 优 位 置 。 

将 式 (11) 代 入 式 (7)， 可 得 最 终 进化 公式 如 式 (12) 所 示 。 


当 5j 越 大 ， 代 
当 取 较 小 值 , 可 防止 粒子 搜索 


Xi t+) =p D+ +0-s,,()) :DD): |C， (一 XDInG 


(12) 


2.3 ”精英 粒子 的 随机 交叉 学 习 

以 上 CE 系数 控制 策略 虽 在 一 定 程度 上 改善 了 QPSO 算法 
搜索 速度 和 收敛 精度 ， 但 随和 迭代 进行 ， 种 群 进化 后 期 多 样 性 不 
可 避免 有 所 损失 ， 使 得 粒子 易 早熟 收敛 。 为 解决 该 问题 ， 本 节 
提出 精英 群体 的 随机 交叉 学 习 策略 ， 使 得 精英 群体 中 粒子 活性 


彭 小 


等 : 基于 


执行 随机 交叉 学 习 ; 


多 根据 式 (12) 更 


改进 QPSO 算法 的 电动 汽车 模糊 控制 器 参 娄 


ChinaXiv 合 作 期 和 


新 每 个 粒子 位 置 ; 


pb) 重复 步骤 pb)~g)， 直 至 满足 一 定 的 循环 结束 条 件 。 
3 ”算法 性 能 分 析 


3.1 


参数 设置 


为 验证 AMF-QPSO 算法 


测试 函数 进行 实验 。 将 AMF-QPSO 算法 分 别 与 固定 系数 、 线性 
递减 、 非 线性 递减 三 类 典 
各 对 比 算法 的 CE 系数 
略 中 xc =0.8; 线性 


低 于 设 定 阔 值 ww 时 , 让 粒子 有 能 力 偏 离 原 位 置 而 进入 其 他 空间 


搜索 ， 使 其 可 能 寻找 到 更 优 位 置 。 其 中 ， 随 机 交叉 学 习 的 定义 
如 下 : 

定义 1 随机 交叉 学 习 。 维 度 为 D 的 搜索 空间 中 ， 令 
和 =00， bh 5) 为 精英 粒子 位 置 ， 为 Dp 维 中 随机 选择 
的 某 一 维 分 量 , 则 X 按 概率 p, 进行 随机 交叉 学 习 后 的 新 位 置 坐 


标 为 X, = (Xs Mop Mp Xpy) » 其 定义 如 下 : 
rand > pp, 


:A 
| 二 ， 其 他 人 


为 保证 算法 收敛 ， p, 采用 自 适 应 线性 递减 策略 ， 如 式 (14) 
所 示 。 

7 (14) 

其 中 : ps 和 pm 分 别 为 最 大 交叉 概率 和 最 小 交叉 概率 ; 1 为 当 


前 代数 ;Gi 为 最 大 和 迭代 次 数 。 
综 上 所 述 , 精英 群体 的 随机 交叉 学 习 伪 代码 见 算法 1 所 示 。 
算法 1 精英 群体 随机 交叉 学 习 


1: ”计算 粒子 适应 度 值 ， 对 粒子 进行 降序 排序 ; 
2: ”选取 适应 度 值 较 优 的 前 16% 个 粒子 组 成 精英 群体 ; 


3: ”计算 精英 群体 中 每 个 粒子 活性 5;, 


om 


a4: if Sij; < Si 

5 按照 交叉 概率 Pw 产生 随机 维 ; 
并 执行 随机 交叉 学 习 ， 如 式 (13); 

6: end if 


2.4 基于 反馈 机 制 的 自 适应 变异 QPSO 算法 
按照 上 述 讨论 ， 反馈 机 制 的 自 适应 变异 QPSO 算法 
(AMF-QPSO) 具 体 执行 流程 如 下 : 
a) 初 始 化 参数 , 包括 种 群 规模 N、 粒 子 旨 
“、 粒 子 搜索 空间 、 初 始 化 空间 等 ; 
b) 在 问题 空间 以 随机 原则 初始 化 粒子 位 置 ; 
9) 按 式 (3) 计 算 粒子 平均 最 优 位 置 ; 


E 度 站 、 最 大 和 代数 


值 上 =0.5 ; 


非 线 1 
Qinia =1.6 ， 12=1.0( 凹 函数 ); 本 文 提 


型 


性 能 ， 利 用 表 2 中 六 个 典型 标 ;# 


CE 系数 控制 策略 进行 对 比分 析 。 


性 递减 CE 


控制 策略 按 如 
递减 CE 系数 控制 策 


固定 系数 控制 策 
各 中 初 值 取 mw =1.0 ， 终 
系数 控制 策略 参考 文献 [9]， 其 中 
出 的 AMF-QPSO 算法 中 ， 


下 设置 : 


经 大 量 实验 验证 ，CE 系数 初 值 mw=08 ，4=0.5 ， 交 又 概率 


Pmax =1.0 ， pm=04 时 ， 算 法 


性 能 最 


佳 
表 2 标准 测试 函数 
围 


函数 名 初始 化 范围 搜索 范围 
Sphere( fi) [=100,50]” [=100,100] 
Rosenbrock( f; ) [-30.15]” [-30,307 
Rastrigin( f ) [-5.12,2.56]” [-5.12,5.12] 
Griewank( f, ) [-600,300]? [=600,600] 
Ackley( f; ) [-32,16]? [-32,32]2 
Schwefel( f; ) [250,500] [-500, 500] 
3.2 结果 分 析 
为 全 面 论 述 算法 性 能 优势 ， 采 用 粒子 规模 N 分 别 为 20、 
50、80， 维 数 刀 分 别 为 10、30、50 维 ， 最 大 迭代 代数 Gmax 分 
别 为 1000、1 500、2 000 代 的 种 群 进行 实验 。 ee 


行 100 次 ， 记 录 结 果 3 


详 见 表 3 所 示 。 


有 反馈 机 种 
象 地 说 明 算 法 性 
Rastrigin 函数 来 观察 


为 


表 3 可 得 ,名 


| 的 


文 提出 的 


适 


色 大 多 数 情况 下 ，CE 系数 
略 所 取得 的 收 和 敛 精 度 均 优 于 其 他 三 种 系数 控 玮 
自 适应 变异 CE 系数 控制 策 
能 优势 ， 
固定 系数 、 
应 变异 四 种 策略 下 适应 度 收敛 曲线 和 CE 系数 的 


并 计算 最 优 解 的 平均 值 。 各 测试 函数 结 


适应 变异 控制 策 
策略 ， 说 明了 具 
略 的 合理 性 与 有 效 性 。 
峰 Sphere 函数 和 多 峰 
线性 递减 、 非 线性 递减 以 及 本 


以 典型 音 


变化 趋势 ， 如 图 2 和 3 所 示 。 


由 图 2 可 得 


现 为 一 个 凸 函 


后 


台 台 
He; 


略 能 根据 粒子 不 同 


期 较 小 3 


和 从， | 


适应 CE 系数 控制 策 
数 , CE 系数 前 


名 在 整个 收敛 过 程 呈 
期 较 大 , 粒子 具有 更 好 全 局 搜索 性 


粒子 具有 更 好 的 局 部 搜索 1 


生 能 ， 自 适应 控制 策 


的 搜索 阶段 来 调整 得 到 


个 合适 的 值 ， 使 得 


算法 具有 更 快 搜 索 速 度 和 更 好 精度 。 此 外 ， 凸 函数 


dj) 计算 粒子 当前 适应 度 值 ， 式 按 (10) 计 算 每 个 粒子 活性 因 
hmm i 
e) 选 择 适 应 度 值 最 优 的 前 I 
(14) 计 算 精 英 粒子 交叉 概率 mm ; 
判断 精英 粒子 活性 , 若 低 于 设 定 下 限 阐 值 s， 


+ 


群体 ， 按 式 


则 按 式 (14) 


同 


万 化 问题 而 不 同 ， 


范围 内 仅 有 


以 这 类 吓 函 


数 


率 较 大 。 


对 于 Sphere 这 类 单 峰 优化 问题 ， 
个 最 优 值 ， 优 化 问题 简单 ， 势 B 


率 针对 不 
因为 全 局 
会 很 快 收敛， 所 


录用 稿 
表 3 不 同 控制 策略 的 实验 结果 
控制 策略 
函数 N D Gmax 
国定 线性 非 线性 自 适 应 
20 10 1000 7.54E-67 3.76E-41 1.64E-27 4.71E-86 
万。 50 30 1500 2.24E-36 3.14E-24 7.68E-20 8.37E-72 
80 50 2000 2.86E-21 1.20E-19 3.97E-13 8.09E-68 
20 10 1000 6.29E+00 5.49E+00 5.28E+00 2.25E+00 
万 50 30 1500 2.27E+01 2.50E+01 2.01E+01 6.41E+00 
80 50 2000 5.35E+01 4.47E+01 4.16E+01 3.87E+01 
20 10 1000 6.60E+00 5.41E+00 9.32E+00 2.98E+00 
f 50 30 1500 6.21E+01 2.08E+01 5.03E+01 1.69E+00 
80 50 2000 1.89E+02 4.47E+01 6.61E+01 1.48E+01 
20 10 1000 9.77E-02 7.79E-02 1.42E-02 9.93E-04 
万 50 30 1500 2.13E-02 1.20E-02 7.05E-03 2.05E-03 
80 50 2000 6.67E-03 5.51E-03 3.35E-03 1.11E-03 
20 10 1000 2.66E-15 2.66E-15 2.45E-15 1.26E-18 
fi 50 30 1500 7.14E-15 5.58E-14 1.19E-10 1.33E-17 
80 50 2000 2.52E-15 3.44E-11 3.25E-08 1.24E-17 
20 10 1000 1.74E+01 1.84E+02 1.68E+00 9.47E-04 
fs 50 30 1500 7.79E+01 3.75E+01 7.83E+00 4.22E-03 
80 50 2000 7.82E+01 1.97E+02 1.45E+01 6.99E-03 
1020 
M1 0 
| 


-100E 一 0— 
10 0 0.2 0.4 0.6 0.8 1.0 ?和 0.2 0.4 0.6 0.8 1.0 
1 x10 1 x10 
(a) 适应 度 曲线 (b) CE 变化 曲线 
图 2 单 峰 Sphere 函数 各 控制 策略 下 适应 度 曲线 


0.2 0.4 0.6 0.8 


1.0 
X 103 


(a) 适应 度 


线 


用 | 


(b) CE 变化 


0.2 0.4 0.6 0.8 


攻 


1.0 
x10 


Da 


局 范 


已 


3 可 得 ， 对 本 


到 3 多 峰 Rastrigin 函数 各 控制 策略 下 适 
F 类似 Rastrigin 的 多 峰 


A 


函数 曲率 相对 较 小 ， 
策略 能 根据 势 凡 


后 


趋势 ， 保 i 


E 粒 子 搜索 前 


期 的 全 局 ] 


内 存在 多 个 极 值 ， 优 化 问题 复杂 ， 势 B 
期 种 群 多 样 性 更 强 。 总 的 来 说 ， 
EH 长 度 相对 变化 趋势 的 不 同 自 
搜索 能 力 和 


j 收 全 


化 问题 ， 由 于 全 


线 
线 


工 


时 
受 ， 


所 以 凸 
自 适应 


适应 调整 CE 变化 


搜索 后 


期 的 局 部 控 


掘 能 


数控 制 策略 为 比较 对 象 ， 与 本 文 提 出 的 具有 反馈 机 


1 


策 


下 


阱 长 度 工 变化 趋势 


自 适 应 变异 控制 策略 进 


省 深层 次 原 


天 


， 以 获得 更 好 收敛 速度 和 搜索 精度 。 


为 探究 本 文 控 仙 ， 以 经 


线性 递减 CE 系 


制 的 CE 系 


行 对 比 ， 观 察 两 种 控制 策 


区 


> 如 


4 所 示 ， 其 


器 


b 


各 下 平均 势 


FE 均 势 阱 长 度 如 式 (15) 


吉 


夜 小 


-25 1 Rastrigin 


A Grie wank 
> A 已 


-35 


0 0.2 0.4 0.6 0.8 1.0 


I 


平 ， 等 : 基于 改进 QPSO 算法 的 


-3 
x10 


(a) 线性 递减 策略 


4 线性 递减 及 


杂 单 峰 优 化 问题 ， 


由 图 4() 可 知 ， 线 性 递减 控制 
单 峰 优化 问题 ， 梯 度 信 


置 ， 种 群 会 迅速 聚 


nr 


9 


函数 


PRastrigin 
A Grie Wank 


50 02 04 06 08 1.0 


1 x10 


(b) 自 适 应 策略 


自 适应 策略 


A 


宋 


日 
势 B 
Il 


明显 ， 
会 短 时 间 


平滑 ， 梯 


息 
4 


搜索 到 全 局 最 优 位 置 ， 


势 阱 收缩 较 组 


个 多 峰 优化 问题 ， 


下 
粒子 会 很 快 搜索 到 
内 收缩 ， Rosenbrock 为 复 


认 
又 


下 势 阱 长 度 变化 
， 类 似 Sphere 的 简单 
全 局 最 优 位 


言 息 不 明显 ， 粒 子 很 难 
曼 ; 而 类 似 Rastrigin 的 四 


度 变 化 情况 均 不 明显 ， 造 成 粒子 在 


局 部 极 值 较 


在 整个 搜索 过 程 势 阱 长 


后 暴 


的 局 部 挖掘 能 力 相 对 较 


弱 ， 粒 子 很 难 搜寻 到 全 局 最 优 位 


EE ， 这 种 控制 


策略 小 算法 收敛 


床 


[ 
单单 峰 


优化 问题 ， 


势 


图 4(b) 可 知 ,在 自 适应 变异 控制 策略 下 ， 类 似 Sphere 的 


速度 ,而且 有 利于 粒子 找 至 
策略 优势 尤为 明显 ， 与 


势 阱 长 度 明显 变 小 ， 势 P 


多 


收缩 更 为 迅速 ， 不 仅 
| 更 优 解 ;类似 Rastrigin 的 多 峰 
4(a) 相 比 ， 
收缩 更 为 剧烈 


9 利于 加 快 收敛 
允 化 ， 
后 期 ， 该 策略 控 
说 明 种 群 进 化 未 


在 搜索 
» 


在 


适应 以 及 精英 粒 了 


4.1 


加 速 


搜索 ， 这 也 是 改进 算法 在 纪 
于 传统 算法 的 原 
基于 改进 QPSO 算 法 的 电动 汽车 调 速 控制 器 设计 
电动 汽车 数学 模型 
蝶 动 汽车 行驶 过 程 中 阻力 主要 来 源 于 滚动 阻力 \ 空 气 阻力 、 
阻力 以 及 坡度 阻力 5 四 


Ba 


o 


个 方面 。 


交叉 学 习作 


色 大 多 数 多 峰 优 


下 ， 粒 子 进 入 小 范 
化 问题 上 收敛 


总 阻力 满足 式 (16) 中 关系 。 


Fa=P, +h,+h, 


其 中 : 深 动 阻力 7 满足 式 (17) 中 关系 。 


F,= fingcosO 


其 中 : nm 为 汽车 


质 


里 ; 


f 为 滚动 摩擦 系数 ; 


数 ， 2 为 行驶 路 面 坡度 。 


空气 阻力 ,满足 式 (18) 中 关系 。 


其 中 : CG, 为 空气 


日 力 系 


+ (16) 


(17) 
8 为 重力 加 速度 系 


(18) 


i 积 ; v 为 汽车 行驶 速度 。 


加 速 阻力 满足 式 (19) 中 关系 。 


其 中 : 5 为 质量 换算 系数 。 
把 坡 阻 力 满足 式 (20) 中 关系 。 


C,Av’ 

“ ”21 坊 

数 ; 4 为 迎风 务 
F =6m 

dt 
F=mgsin0 


(19) 


(20) 
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汽车 行驶 过 程 中 ， 行 驶 速度 与 车 轮转 速 之 间 满 足 式 (21) 关 
系 。 


v=0377- CD 
其 中 : 为 电机 转速 ，7 为 车 轮 半 径 ， i 为 变速 器 传动 比 : 为 
主 减速 器 传动 比 。 
汽车 行驶 过 程 中 ， 总 转 矩 和 所 受 总 阻力 满足 式 (22) 关 系 。 
Ts = Fer (22) 
有 机 负载 转 矩 与 总 转 失 Ts 之 间 满 足 
ee 
i (23) 
其 中 : 7 为 传动 效率 。 


4.2 电机 数学 模型 
动 汽车 以 BLDCM 为 驱动 器 051， 为 建立 电机 数学 模型 ， 
进行 如 下 假设 : 

a) 不 考虑 磁 清 损耗 与 涡流 损耗 ; 

b) 忽略 磁 路 饱和 ， 和 忽略 磁 横 效应 ; 

c) 三 相 绕组 均匀 对 称 分 布 于 电 枢 绕组 表面 。 
基 尔 霍 夫 电压 定律 可 得 ， 任 意 一 相 绕组 上 电压 平衡 方程 
如 式 (24) 所 示 。 


口 = 应 +L 策 + 下 C4) 
dt 


对 于 三 相 绕 组 ， 可 等 效 如 式 (25) 所 示 。 


ul IR 0 Of Eh Ly Fel i le 
wl=I0 ER Obltib, Lb, blpli |+tle, (2 5 ) 
U, 0 0 Ri Ls L, LL i € 


其 中 : wu、w、u 为 三 相 定 子 电 压 (V); 6、4、e 为 三 相 定 子 反 电 
动 势 (V); i 认为 三 相 定子 相 电 流 (A); 六 乒 到 为 三 相 定子 
自 感 ( 外， 局 bb bh、 LL 为 三 相 定子 绕组 之 间 的 互感 
(四 ; RR、B、R, 为 三 相 定子 绕组 的 相 电 阻 (Q); 2 为 微分 算 子 
(dih。 


为 便于 分 析 ， 现 假设 三 相 绕组 完全 对 称 分 布 ， 且 绕组 连接 
方式 采用 星 型 无 中 线 式 ， 则 有 
五 = 五 = 到 = 
R,=R,=R=R 
(26) 


Lys=L. = = =L,=L,=M 
iiti,t+i.=0 


其 中 : 工 和 AM 分 别 为 定子 自 感 和 互感 ， 则 式 (25) 可 写成 


有 RO 0 L-M 0 0 人 e。 
u,|=I0 R 0 +| 0 L-M 0 |pli,l+le,| (27) 
u, 0 0 RR 0 0 L-M i. e. 


罗 


I 
c 


式 (27) 可 得 BLDCM 的 等 效 电 路 如 图 5 所 示 。 
fa L-M 
Ua VY YY 
lp L-M 
Up 
ER 


图 5 BLDCM 等 效 电 路 图 
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运行 状态 下 ，BLDCM 的 感应 电动 势 与 转速 、 绕 组 古 数 成 
正比 ， 因 此 电 枢 绕组 反 向 电动 势 的 公式 可 用 式 (28) 表 示 。 


三 

E= Tb =K, (28) 
其 中 : 2p, 为 极 对 数 ，N 表示 各 相 线圈 串联 臣 数 ;办 表示 每 极 磁 
通 量 ; a 为 极 弧 系数 ; 7 为 电机 转速 ; K, 为 反 电 动 势 系数 


(Krpm/V )。 
有 机 的 电磁 转 矩 方程 和 运动 方程 如 式 (29) 所 示 。 


T= (e +e,i, +e.i,) 
ee 


n 
更 (29) 
J—=T-T.-Bn 
dt 


其 中 : 对 于 单 相 ,可 得 =Ki*i, ; J 为 转子 惯量 (gcem”); 8 为 阻 
尼 系 数 ，Ki 为 转 矩 常数 ( N.m/4); 7 为 电机 转速 ， 为 电机 负 
载 转 矩 ( Nm)， 由 式 (23) 决 定 。 

对 上 述 方 程 组 取 拉 普 拉 斯 变换 , 可 得 BLDCM 的 动态 模型 ， 
如 图 6 所 示 。 


7,0) 
U,(s) 天” 二 > 1 2s) 
Ls+R Js 
加 
图 6 BLDCM 动态 模型 图 
由 图 6 可 得 ， 电 机 的 传递 函数 可 用 式 (30) 表 示 。 
K7 LS+R 

MIT RS RR J RS KR) (G30) 


4.3 ”基于 AMF-QPSO 算法 的 控制 器 设计 
动 汽车 行驶 工 况 复杂 ， 且 需 进行 频繁 调 速 ， 传 统 PID 控 
制 器 很 难 取得 预期 控制 效果 。 本 文采 用 模糊 PID 作为 控制 器 ng 
系统 结构 如 图 7 所 示 。 


速度 检测 速度 检测 


给 定 


模糊 PID 
控制 器 


电流 检测 RR 一 一 


图 7 基于 模糊 自 适 应 PID 的 BLDCM 调 速 系统 
影响 控制 器 性 能 指标 的 主要 有 两 大 因素 : a) PID 控制 器 参 
数 &、k、ks; b) 量化 因子 上 、K 和 比例 因子 如。 由 于 模糊 
PID 控制 器 对 于 k,。、kh 、ko 具有 一 定 调节 能 力 ， 但 控制 器 本 身 
自 适应 性 能 不 高 ， 为 进一步 优化 模糊 控制 器 性 能 ， 可 从 模糊 控 

制 器 关键 参数 ， 即 量化 因子 和 比例 因子 的 优化 入 手 。 
若 以 被 控 量 误差 及 其 误差 变化 率 ec 作为 系统 输入 ， 设 。 
的 基本 论 域 为 [ei,en]，ec 的 基本 论 域 为 [eci,ecn], 输出 变量 4 的 
基本 论 域 为 [wm 。 误 差 变 量 模 糊 论 域 
误差 变量 率 的 模糊 论 域 


-jtl…0wj-lj,jeN ， 输 出 控制 量 模糊 论 域 


ER 


iths05ili},ieN ， 


吊 


{ 
{ 


谤 
和 
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{kk+b…,0,…,k 一 b 旭 ,keN， 则 量化 因子 、K. 和 比例 因子 态 重 的 性 能 指标 来 确定 大 小 。 
式 (31) 确 定 。 
5 ”实验 与 结果 分 析 
2i 2j Up 一 三 
k k= 起 尖 31 
er 一 er ecr 一 eCr 2k ( ) 5.1 搭建 实验 平台 

但 环境 因素 及 各 项 参数 误差 ， 使 得 理论 值 与 实际 值 常 存在 实际 电动 汽车 调 速 控制 系统 工 况 复杂 ,包括 软 、 硬 件 系 统 。 
偏差 ， 造 成 系统 性 能 随机 性 强 ， 对 控制 参数 敏感 ， 所 以 实际 应 。” ”为 便于 实验 室 研究 , 硬件 系统 设计 上 ， 本 文 以 Maxon 公司 制造 
用 时 还 需 根据 经 验 进行 调整 。 的 小 型 BLDCM 模拟 电动 汽车 电机 。 查 阅 数据 手册 ， 详 细 参 数 

因此 ,本 节 以 此 为 研究 重点 ,在 模糊 PID 基础 上 ,利用 AMF-  ” 见 表 5 所 示 。 以 ST 公司 的 专用 集成 芯片 L6235 为 核心 搭建 逆 


QPSO 算法 优化 模糊 PID 控制 器 中 量化 因子 、K 和 比例 因子 。 变 电 路 。 


的 整定 过 程 ， 使 其 能 适应 系统 参数 变化 ， 从 而 增强 系统 抗 干 表 5 直流 无 刷 电机 参数 表 
扰 和 自 适 应 性 能 。 算 法 结构 如 图 8 所 示 。 主要 指标 参数 值 
/AMF-QPSO 算 法 (Simulink 模型 额定 电压 /U 24V 
人 额定 功率 /P 90 W 
一 一 >| 初始 化 粒子 粒子 依次 赋值 给 额定 转速 人 h 5700 r/min 
粒子 群 | Ks ks hk 额定 转 矩 /T 47mN:m 
人 i 相 感 抗 人 1.6x105H 
TIN 相 电阻 /有 R 0.51Q 
生生 满 证 名 区 反 电 动 势 /KE 1.47x103 V.s/rad 
止 条 件 Lem ee 转 矩 常数 /Kr 0.014 Nma/A 
转子 惯量 / 5.54x107 Kg-m? 
~ 7、 J 软件 设计 上 ， 利 用 ActiveX 技术 将 LabVIEW 和 MATLAB 
图 8 AMF-QPSO 算法 优化 模糊 PID 联系 起 来 。LabVIEW 主要 完成 数据 采集 与 处 理 、 系 统 拟 合 
整个 系统 可 分 为 AMF-QPSO 算法 执行 单元 和 Simulink 仿 ” PWM 控制 信号 输出 、 网 络 通信 以 及 图 形 化 界面 设计 [3]， 


单元 两 大 部 分 。 其 中 AMF-QPSO 算法 由 MATLAB 中 的 S- MATLAB 用 于 执行 AMF-QPSO 算法 以 及 对 比 算法 。 整 个 实验 
unction 函数 实现 。 设 粒子 维 数 为 3， 分 别 代表 2 个 量化 因子 ”平台 系统 框架 如 图 9 所 示 。 
和 1 个 输出 因子 ， 种 群 规模 为 20， 最 大 迭代 代数 为 500， 其 他 
控制 参数 参考 3.1 小 节 。 模 糊 PID 控制 器 及 BLDCM 控制 系统 二 == 和 
Simulink 工具 箱 搭建 , 其 中 模糊 PID 控制 器 以 BLDCM 控制 | 入 
系统 转速 误差 。 及 其 变化 率 ec 作为 模糊 控制 器 输入 ， 以 调整 量 ER 二 二 ee 


es 


Ak。、Ak, 、Ako 为 控制 器 输出 ,模糊 集 为 {NB、NM、NS、ZE、 
PS、PM、PB}。 各 变量 基本 论 域 和 模糊 论 域 如 表 4 所 示 。 模 糊 9 实验 平台 系统 框图 
规则 与 文献 [17] 相 同 ， 丸 属 度 函 数 采用 三 角 型 ， 以 if…… then ”5.2 控制 器 动态 性 能 分 析 
的 形式 输入 模糊 控制 规则 ， 运 用 重心 法 将 模糊 量 进行 清晰 化 处 为 验证 控制 器 动态 性 能 ,采用 单位 阶 跃 信号 作为 测试 信号 ， 
里 得 到 被 控 系 统 控 制 量 。 MATLAB 平台 上 依次 运行 传统 模糊 PID(Fuzzy) 算 法 、QPSO 算 
表 4 模糊 自 适应 PID 控制 器 参数 表 法 优化 后 的 模糊 PID(QPSO-Fuzzy) 算 法 以 及 AMF-QPSO 算法 
参数 e ec Ak; Ak, Ab 优化 后 的 模糊 PID(AMF-QPSO-Fuzzy) 算 法 。LabVIEW 利用 
语言 变量 E EC ks k k, ActiveX 技术 调用 MATLAB 算法 进程 , 并 通过 数据 采集 卡 (DAQ) 
基本 论 域 [-100 100] [-105109] [-0.01 0.01] [-0.1 0.1] [-0.01 0.01] 获取 的 电机 实际 工作 状态 ， 根 据 4.1 和 4.2 节 构 建 的 数学 模型 
模糊 论 域 。 [-6 6] [-6 6] [-6 6] [-6 6] [-6 6] 拟 合 电 动 汽 车 速度 控制 系统 模型 ， 最 后 输出 PWM 控制 信号 ， 
转换 因子 。 0.06 6/105 1/600 1/60 1/600 控制 道 变 电路 。 在 上 位 机 界面 上 观察 系统 输出 响应 曲线 ， 实 验 
为 平衡 系统 动 .静态 特性 , 在 保证 稳 态 误差 趋 于 零 的 同时 ， 结果 如 图 10 所 示 。 
使 系统 的 超 调 量 o、 调 节 时 间 4 控制 在 较 小 范围 内 ， 则 适应 度 为 定量 比较 优化 后 的 控制 器 的 性 能 优势 ， 计 算 各 控制 系统 
函数 确定 为 单位 阶 跃 下 的 性 能 指标 ,包括 稳 态 误差 、 超 调 量 以 及 调节 时 间 ， 
fimess = | kle(Ol dr + ht, + ho (32) ee | 
: 根据 图 10 及 表 6 实验 数据 可 知 ， 经 QPSO 及 AMF-QPSO 
其 中 : 太 、 已 、 有 为 各 性 能 指标 的 权重 系数 ， 由 实际 系统 所 侧 ” 算法 优化 后 的 控制 器 在 稳 态 误差 、 超 调 量 以 及 调节 时 间 上 均 优 
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于 传统 模糊 PID 控制 器 。 经 AMF-QPSO 算法 优化 后 的 控制 器 优势。 


较 QPSO 算法 优化 后 的 控制 器 复杂 程度 更 高 ， 故 在 上 升 时 间 上 ”5.3 ”控制 器 调 速 性 能 分 析 

略 慢 与 QPSO 算法 优化 后 的 控制 器 ， 但 AMF-QPSO 算法 优化 对 系统 调 速 性 能 进行 测试 ， 仿 真 时 长 0.4s， 主 要 通过 空 载 

后 的 控制 器 在 稳 态 误差 、 超 调 量 以 及 调节 时 间 上 具有 绝对 优势 。 “和 负载 实验 进行 论述 。 系 统 空 载 状况 下 启动 ， 设 定 预期 转速 为 
1000t/min, 速度 稳定 后 , 在 t=0.1s 时 , 突 加 负载 7=0.03N:m， 


Ei 模拟 行驶 坡度 变 为 5”， 坡 度 阻力 相应 增加 时 系统 速度 响应 ; 
= | 在 {=-0.2s 时 ， 将 转速 调节 至 2 000 rmin， 模 拟 电动 汽车 突然 加 
| 速 时 系统 速度 响应 ， 观 察 系 统 速 度 响应 曲线 、 转 速 及 其 局 部 放 
大 响应 曲线 如 图 13 和 14 所 示 ， 转 盾 响 应 曲线 如 图 15 所 示 。 


转速 (n/min) 


图 10 各 控制 策略 下 的 阶 跃 响应 曲线 | 

表 6 各 控制 系统 单位 阶 跃 下 的 性 能 指标 | L 

< 二 优化 算法 稳 态 误 差 5/%。” 超 调 0/% 调节 时 间 4/s I 

0 0.08 0.16 0.24 0.32 0.4 
Fuzzy 0.22 18.5 0.036 t/s 
QPSO-Fuzzy 0.18 8.7 0.022 图 13 各 控制 策略 下 的 速度 响应 曲线 

AMF-QPSO-Fuzzy 0.08 0.5 0.012 由 图 13 可 知 ，AMF-QPSO 算法 优化 后 的 控制 器 在 调 速 性 
为 进一步 验证 AMF-QPSO 算法 优化 后 的 模糊 PID 控制 器 ”能 上 优 于 传统 模糊 PID 控制 器 及 QPSO 算法 优化 后 的 模糊 PID 


的 抗 干扰 及 自 适应 性 能 ， 外 加 幅 值 为 1 的 方 波 信号 ， 模 拟 电动 控制 器 ， 其 性 能 优势 详 见 图 14 所 示 。 
车 行驶 过 程 中 受到 的 周期 干扰 信号 。 实 验 结果 如 图 11 所 示 ， 
节 部 分 如 图 12 月 


1 


0.0196 0.02 


be 
i -0.5 
1004 
1.5 ye QPSO Fyzzy ! 
全 1002 
写 1 
四 | 
浸 1000 
' 
1.0158 | 
; 998 上 
| 0.0999 0.1 0.1001 0.1002 0.1003 0.1004 
1 t/s 
1 
> 0.7158 上 -一 -一 yy 多 ------ | (b) 突 加 负载 时 速度 响应 
一 站 一 AMF-QPSQ-Fuzzy | 
—et— Fuzzy 2010 Be 1 1 
QPSO-Fyzzy 一 x | QPSO-Huzzy 1 | 
0.4158 广 003 | 
在 | 
0.4 0.6 0.8 六 2000 
t/s 
1995 
图 12 各 控制 策略 下 的 方 波 响 应 曲线 局 部 图 
1990' v 
图 11、12 可 知 ，AMF-QPSO 算法 优化 后 的 模糊 PID 控 Wa Om 0 90 O39 
制 器 虽 在 响应 速度 上 略 慢 与 QPSO 算法 优化 后 的 控制 器 ， 但 控 (c) 稳 态 时 速度 响应 
制 器 具有 超 调 小 、 稳 态 精 度 高 、 抗 干扰 及 自 适 应 性 能 强 等 绝对 图 14 各 控制 策略 下 速度 响应 局 部 曲线 
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到 


14(a) 为 各 控制 系统 空 载 启 动 时 速度 响应 ， 可 得 经 AMF- ”加 快 后 期 的 搜索 速度 和 收敛 精度 。 为 验证 算法 性 能 ， 将 AMF- 

QPSO 优化 后 的 控制 器 在 响应 速度 上 虽 略 慢 于 QPSO 优化 后 的 ”QPSO 算法 应 用 于 电动 汽车 调 速 控制 系统 中 模糊 控制 器 参数 优 

空 制 器 ， 但 在 稳 态 精度 上 最 优 ， 误 差 仅 为 0.05%;”(b) 为 t=0.1s ”化 问题 中 ， 并 以 LabVIEW 为 实验 平台 ， 通 过 与 传统 模糊 PID 

时 ， 突 加 负载 后 各 系统 速度 响应 ， 用 于 验证 系统 抗 干 扰 能 力 。 控制 器 以 及 QPSO 优化 后 的 模糊 PID 控制 器 进行 性 能 对 比 ， 实 

里 论 知识 可 知 ， 突 加 负载 后 电机 转速 会 短 时 下 降 ， 此 时 控制 仿 结 果 表 明 AMF-QPSO 算法 优化 后 的 控制 器 具有 良好 的 响应 

器 为 维持 速度 恒定 ， 会 迅速 提升 绕组 电流 ， 进 而 增 大 输出 转身 速度 和 控制 精度 ， 且 对 于 系统 扰动 具有 很 好 的 自 适 应 和 自 调节 

直至 平衡 负载 转 矩 ， 使 得 系统 能 尽快 跟随 系统 参数 变化 。 其 中 能 力 ， 对 电动 汽车 速度 控制 器 优化 具有 重要 意义 。 
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